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sinvn^RY 

Modern  aircraft  design  trends  have  been  toward 
aircraft  with  long,  dense  fuselages,  smaller  aero- 
dynamic surfaces,  swept  wings  and  higher  wing  loadings » 
These  aircraft  are  required  to  fly  at  high  altitudeso 
Flight  data  from  such  aircraft  have  indicated  that  the 
linearized  dynamic  equations  do  not  always  predict  the 
dynamic  behavior  accurately.   In  addition,  certain  cross- 
coupling  effects  were  developing  between  the  lateral 
and  longitudinal  motion  of  the  aircraft. 

The  effect  of  the  trend  toward  concentrating  the 
mass  in  the  fuselage  was  investigated  by  making  an 
engineering  approximation  of  the  change  in  the  moments 
of  inertia  caused  by  this  change  in  mass  distributiono 
Four  configurations  were  considered,  each  one  approach-" 
ing  nearer  the  limit  of  a  missile  or  body  of  revolutiono 
Since  the  lower  lift  curve  slopes  and  smaller  wing  areas 
of  the  swept-wing  aircraft  sometimes  require  flight  at 
very  high  angles  of  attack,  the  effect  of  inclination 
of  the  principal  longitudinal  axis  to  the  flight  path 
was  investigated.   The  effect  of  sweep  was  approximated 
by  varying  the  effective  dihedral  parameter,   C.        „ 
In  addition  the  effect  of  varying   C„   ,  the  "weather- 
cock-stability"  parameter,  was  studied  since  long  fuse- 
lages affect  this  parameter  adversely. 

The  investigation  was  made  using  a  Reeves  Electronic 
Analogue  Computer  (REAC)  by  setting  up  a  mathematical 
model  of  a  representative  military-type  airplane  flying 
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at  a  Mach  number  of  0„6  at  sea  level  and  at  40,000 
feeto   The  mathematical  model  was  based  on  equations 
of  motion  v\rhich  were  extended  to  include  products  of 
the  perturbation  velocity  termso   The  effect  of  changing 
the  mass  distribution,  the  inclination  of  the  longi- 
tudinal principal  axis,  and  variation  of  C^   and 
Ci        was  investigated  by  changing  the  mathematical 
model  as  required  by  the  parameter  change  considered 
and  obtaining  the  REAC  response  to  a  fixed  or  standard 
step-function  control  input* 

The  results  indicated  that  the  lateral  motion 
coupled  into  the  longitudinal  motion  and  that  this 
effect  was  increased  as  the  fuselage  became  longer 
and  denser.  Altitude  aggravated  this  coupling  as  it 
did  any  trend  toward  instability  caused  by  a  change  in 
the  other  parameters  consideredo   The  rudder  caused 
greater  coupling  than  the  ailerono   A  positive  inclina- 
tion of  the  principal  longitudinal  axis  caused  an 
increase  in  dynamic  stability  and  a  decrease  in  coupling 
viiile  a  negative  inclination  had  the  opposite  effecto 
A  decrease  in  Q^  produced  a  marked  trend  toward  in- 
stability  and  an  increase  in  coupling  while  a  decrease 
in   G   had  the  opposite  effecto 

The  investigation  required  by  this  thesis  was 
carried  out  in  the  Computer  Center  of  the  University 
of  Minnesota* 
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INTRODUCTION 

Fighter  aircraft  of  recent  design  show  a  trend 
toward  long  fuselages  where  most  of  the  mass  is  con- 
centrated, smaller  aerodynamic  surfaces,  high  rates  of 
roll,  and  increased  use  of  sweepback.   Flight  data 
from  such  aircraft  indicate  that  the  linear  equations 
of  motion,  which  result  after  several  simplifying 
assumptions,  do  not  always  predict  accurately  the 
dynamic  behavioro   In  addition,  these  same  design 
trends  have  required  more  emphasis  to  be  placed  on  the 
consideration  of  the  dynamic  stability  of  an  aircraft 
in  the  design  stageo   In  earlier  practice,  dynamic 
stability  usually  took  care  of  itself  after  the  static 
stability  requirements  were  met. 

The  effect  on  coupling  and  dynamic  stability 
caused  by  the  parameters  principally  affected  by  these 
design  trends  was  investigated  by  means  of  a  mathe- 
matical model  entered  on  the  Reeves  Electronic  Analogue 
Computer  (REAC)o   An  attempt  was  made  to  use  a  repreaenta- 
tive  or  typical  aircraft  of  the  fighter  type*  The 
mathematical  model  was  based  on  equations  of  motion 
vrtiich  were  extended  to  include  the  products  of  the 
velocity  change  usually  assumed  to  be  zeroo 

The  mass  distribution  change  was  investigated  by 
considering  four  configurations  each  more  nearly 
approaching  a  missile  from  the  standpoint  of  its  moment 
of  inertia  propertieSo   The  effect  of  the  inclination 
of  the  principal  longitudinal  axes  was  considered  at 
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inclination  angles  of  zero,  plus  ten,  and  minus  ten 
degreeSo   The  effect  of  Q   and  Ca      was  studied  by 
changing  the  mathematical  model  as  required  by  a  change 
in  these  parameters  over  a  range  of  values.   Altitude 
effects  were  considered  by  calculating  the  dimensional 
derivatives  for  40,000  feet,  changing  the  mathematical 
model  accordingly,  and  running  the  program  as  it  was 
done  for  the  sea  level  case. 

The  effect  of  a  parameter  change  was  observed  by 
obtaining  successive  REAC  responses  to  a  standard  step 
function  control  input*   Only  the  parameter  investigated 
was  changed  so  that  any  variation  in  the  time  history 
response  was  attributed  to  the  parameter  changeo   The 
data  was  presented  by  superimposing  these  successive 
time  historieSo   The  trend  of  the  effect  of  the  para- 
meter change  considered  could  then  be  easily  observedo 
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DEVELOPMENT  OF  THE  PROBLEM 
The  classical  approach  to  airplane  dynamics  is 
outlined  here  in  a  very  brief  raannero   The  airplane 
was  considered  to  be  a  rigid  body«   Six  equations 
were  required  to  define  its  motion;  one  force  equation 
and  one  moment  equation  for  each  axis«   These  equations 
included  the  force  and  moment  components  due  to  aero- 
dynamic and  mass  reactions  together  with  the  force  of 
gravity  and  thrust  effects o   Second  order  effects  were 
neglected  and  the  equations  simplified  to  a  degree 
which  permitted  hand  solution. 

Appendix  I  develops  the  equations  of  motion 
used  in  this  investigatiouo   The  classical  approach 
has  been  used  in  general  but  certain  departures  have 
been  made*,   In  Appendix  II  the  equations  developed  in 
Appendix  I  have  been  simplified  by  assuming  that  the 
airplane  initially  is  in  straight  and  level,  steady, 
gliding  flight.   The  thrust  effects  were  allowed  to  go 
to  zero  since  they  are  not  critically  essential  in  a 
dynamic  stability  studyo 

Derivative  terms  of  order  higher  than  one  in  the 
series  expansion  of  the  expressions  for  the  forces  and 
moments  were  neglectedo   The  effect  of  these  terms  is 
known  to  be  secondary,  data  was  available  only  for  the 
first  derivative,  and  inclusion  of  these  terms,  even  if 
available,  would  require  more  components  than  the  REACTS 
provided,   A  second  assumption  generally  made  is  that 
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the  product  of  the  perturbation  velocities  is  essenti- 
ally zeroo   This  yields  two  independent  sets  of  three 
linear  equationso   This  assumption  was  not  madeo   The 
equations  simplified  in  Appendix  II  are  presented  belowt 


■^  -p  lift   +-  -|?a  f  ^tj-^x  ]  -  <fR.  Dv 
^  ^t  /        li 

The  axes  used  and  the  sign  convention  observed  is 
presented  in  Table  I  and  illustrated  in  Fig.  S  ^(Aiich 
heads  Table  I.« 

Refo  2  points  out  that  motion  of  the  longitudinal 
symmetric  type  does  not  induce  motion  of  the  lateral 
assymmetric  type.   An  examination  of  the  equation  set 
presented  above  indicates  that  there  are  no  terms  vdiich 
would  couple  symmetric  motion  into  the  assymmetric 
motion  if  the  assymmetric  motion  were  initially  zerOo 

Duncan  in  Refo  2  also  points  out  that  when  devia= 
tions  are  so  small  that  their  squares  and  products  can 
be  neglected,  there  will  be  no  coupling  of  the  lateral 
motion  into  the  symmetric  motion©  When  the  deviations 
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are  not   small   longitudinal  motion  will  persist   alone, 
but   a  purely  lateral  motion  will,    in  general,    induce 
longitudinal  motiono      The   motion  will   couple   into  the 
forward  velocity  through  the  term   'lt/l  5   the   pitching 
velocity  through  the  terms    1?ri.  f  -  x  - 1 1- ]    n^   J%z  and 

"p"^  !><*-       ;    and  into  the   normal   velocity,    through  the 

h 

term  ipr    o   It  can  be  observed  that  the  coupling  terms 
are  modified  generally  by  moments  and  products  of  inertia, 

The  study  of  the  effect  of  the  trend  in  mass  dis- 
tribution toward  a  long,  dense  fuselage  was  made  by 
considering  four  cases  starting  with  a  hypothetical 
fighter  airplane  of  conventional  typeo   The  character- 
istics of  this  basic  airplane  are  presented  in  Table  IIo 
The  fuselage  was  considered  to  represent  60  per  cent 
of  the  airplane *s  mass.   It  was  taken  to  have  a  maximum 
width  of  five  feet  and  to  have  an  eight-to°one  slender- 
ness  ratiOo  The  pitching  moment  of  inertia  calculated 
using  these  assumptions  and  the  approximation  that,  for 
purposes  of  making  an  engineering  estimate,  the  fuse- 
lage can  be  considered  a  long,  slender  cylinder  resulted 
in  a  calculated  value  of  the  pitching  moment  of  inertia 
which  was  very  near  the  Ii^  of  the  airplane.   As  the 
mass  became  completely  concentrated  in  the  fuselage  it 
was  assumed  that  the  slenderness  ratio  would  increase 
to  twelve  to  one,  the  fuselage  width  would  remain  the 
same,  and  that  the  pitching  and  yawing  moments  of 
inertia  could  be  calculated  using  the  long  cylinder 
approxiraation.   In  this  final  form,  the  aircraft  would 


be  a  missile  with  1  g.  equal  to  Tc|  and  X^  very  small, 

Figo  10  illustrates  the  engineering  approximation 
that  was  made  of  the  variation  in  moments  of  inertia 
as  the  fuselage  became  longer  and  denser.  The  mass 
distribution  changes  were  made  v\^ile  the  aerodynamic 
damping  terms,  control  power  terms,  mass  and  other 
effects  were  considered  constant.   This  would,  of 
course,  be  unrealistic  but  it  has  the  advantage  of 
studying  the  effect  of  a  change  in  one  parameter  in 
the  "partial  derivative"  sense. 

The  product  of  inertia  term  resulted  largely  from 
the  fact  that  stability  axes  were  used.  When  the 

•)C  -axis  is  so  oriented  in  an  airplane  in  the  steady 
flight  condition  that  it  is  parallel  to  the  relative 
wind,  the  Eulerian  axes  are  referred  to  as  stability 
axes,  NACA  TR  5^9  indicated  that  the  T^^   terms  were 
unimportant  for  conventional  aircraft,  however,  Ref.  7 
showed  that  the  I^a  term  may  have  a  pronounced  effect 
on  the  lateral  stability  of  some  high  speed  airplanes. 
It  has  already  been  noted  that  the  Ty^  term  is  involved 
in  the  coupling  of  the  lateral  and  longitudinal  motion. 
If  the  longitudinal  axes  are  inclined  at  an  angle 
Y[  to  the  flight  path,  the  moments  of  inertia  about  the 
flight  path  axis  and  the  axis  normal  to  the  flight 
path  are: 

Ix.  '   'Clio"  ^^J  5^^^5)^ 
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The  angle  \]  largely  results  from  the  wing  incidence 
angle  and  the  angle  of  attacko 

The  trend  toward  swept  wings  and  long  fuselages 
makes  C^  and  C^  increasingly  iraportanto   Long  fuse- 
lages adversely  affect  C^  ,  the  "weather-cock" 
parameter o   Swept  wings  have  very  large  values  of 
effective  dihedral  thus  affecting,   C^^   ,  the 
effective  dihedral  parametero   The  effect  of  a  change 
in  these  parameters  on  coupling  and  on  dynamic  stability 
was  investigatedo  The  change  in  the  static  derivative 
characteristics  was  selected  for  study  since  general 
flying  characteristics  are  influenced  primarily  by  the 
static  stability  parameters  according  to  Refo  60 

The  coefficients  such  as  L ^  are  dimensional 
stability  derivativeSo   They  are  direct  functions  of 
airframe  geometry,  mass,  mass  distribution  and  flight 
conditions;  in  addition  to  being  direct  functions  of 
the  non-dimensional  stability  derivatives  such  as, 
Qn       „  These  basic  non-dimensional  stability  derivatives 
are,  in  turn,  implicit  functions  of  airframe  geometry, 
mass,  mass  distribution  and  flight  conditionso   The 
values  of  the  non-dimensional  stability  derivatives 
for  the  basic  airplane  appear  as  Table  II o  The  dimen- 
sional stability  derivatives  are  defined  in  Tables  III 
and  IVo   The  values  of  the  dimensional  derivatives  were 
calculated  and  appear  as  Table  Vo  RefSo  1,  2,  3?  k^    7 9 
8,  and  9  discuss  these  derivatives,  their  characteristics, 
and  how  they  are  determinedo   They  are  assumed  to  be 
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constant  for  a  certain  range  of  applicability  and  flight 
conditionSo 

The  dimensional  stability  derivatives  can  be  used 
directly  as  numerical  coefficients  in  the  dynamic 
equations  of  the  airframe  based  on  real  timeo  Thus, 
stability  derivatives  in  this  form  are  used  in  setting 
up  the  mathematical  model  for  the  analogue  computero 
The  use  of  this  form  of  derivative  permits  the  exten- 
sion of  the  non-dimensional  derivatives  to  other  flight 
conditions  vd.th  only  secondary  error o  This  was  done 
to  obtain  the  mathematical  model  of  the  airplane  at 
40,000  feeto   Certain  values  of  the  non-dimensional 
derivatives  \f\^ich  were  specifically  applicable  to 
flight  at  this  altitude  were  used  while  others  were 
extended  by  the  change  in  the  values  of  p        and   U   o 

Phillips  in  Refo  4  considered  the  effect  of 
steady  rolling  on  longitudinal  and  directional  stabil- 
ityo  The  equations  v*iich  he  used  can  be  arrived  at 
from  the  equations  developed  in  Appendix  I  by  applying 
the  assumptions  that  he  made  in  order  to  arrive  at  an 
equation  set  which  could  be  handled  analy1:icallyo 
He  retained  the  product  terms  which  couple  the  longi- 
tudinal and  lateral  mot  ions »   He  found  that  rapid 
rates  of  roll  caused  instability  if  the  directional 
and  longitudinal  stabilities  are  different  when  the 
rolling  frequency  exceeds  the  lower  of  the  pitching 
and  yawing  natural  frequencies  of  the  non-rolling 
airplaneo   The  instability  existed  only  v^ile  the  air- 
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plane  was  rollingo   In  airplanes  of  short  span  and 
high  fuselage  density  flying  at  high  altitudes, 
this  instability  may  cause  dangerous  attitude 
changes o   Further,  he  found  that  vhen  an  airplane 
rolls  about  an  axis  vvhi  ch  is  not  alligned  with  its 
longitudinal  principal  axis,  forces  were  introduced 
which  tended  to  swing  the  fuselage  out  of  line  with 
the  flight  patho   These  are  the  inertial  forces 
usually  neglected  Wien  the  linearized  equations  are 
used. 

In  Refo  4,  Phillips  made  many  assumptions  in 
an  effort  to  obtain  equations  >hi  ch  were  capable  of 
hand  solution  yet  would  predict  the  effect  on  coupling 
of  steady  rolling.   In  the  investigation  conducted 
incident  to  this  thesis,  the  effects  of  the  three  motions, 
roll,  pitch,  and  yaw,  on  coupling  as  produced  by  step^ 
function  control  inputs  were  investigatedo   As  many 
influencing  factors  were  retained  as  the  REAC*s  compon- 
ents would  permit* 

Solutions  of  the  equations  of  motion  used  in  this 
thesis  would  require  the  simultaneous  solution  of  six 
differential  equations,  three  of  second  order,  made 
non-linear  by  product  terms  such  as  an  and  -pre     o 
The  labor  involved  would  have  been  prohibitive*  The 
equations  were  fitted  to  two  REACTS  and  four  servo- 
raultiplierSo   The  solutions  were  recorded  simultane- 
ously on  two  four-channel  Sanborn  recorders© 


-12- 

EQUIPMENT  AND  PROCEDURE 

The  equipment  used  to  conduct  the  investigation 
consisted  of  two  REAC's,  four  servo-multiplier  units, 
two  four  channel  Sanborn  recorders,  and  a  one- shot 
multi-vibrator  time  switch*   The  entire  computing 
system  is  illustrated  in  the  photograph  presented  as 
Fig.  2o   The  Sanborn  recorder  is  shown  in  Fig o  3,  the 
lateral  motion  and  control  function  computer  in  Figo  4j 
the  longitudinal  motion  computer  in  Figo  5,  and  the 
time  switch  in  Fig«  6o 

Each  of  the  four  servo-mechanisms  drove  three 
multiplying  potentiometers o   These  multipliers  provided 
(^Q   of  the  product  of  the  inputs o   This  restriction 
made  it  necessary  to  build  up  the  perturbation  velocities 
to  ten  times  their  actual  value  in  order  to  obtain  a 
simple  product  such  as  -pfi,  o      The  additional  amplifi- 
cation required  by  the  servo-multipliers  made  computer 
balance  very  delicate  and  this  was  a  source  of  difficulty. 

The  Sanborn  recorder  made  time  history  records  of 
four  quantities  simultaneously.   The  imprint  was  made 
by  a  heated  pen  on  sensitized  papero   The  point  did  not 
touch  the  paper,  thus  eliminating  the  drag  usually 
associated  with  a  recording  pen. 

The  time  switch  was  designed  to  provide  an  accurately 
timed  step-function  as  a  standard  input »   Basically,  the 
switch  is  a  one-shot  multi-vibrator  Wiich  fires  just 
once  viien  the  operating  button  is  pressed*   This  allows 
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current  to  flow  in  the  coil  of  a  relay  in  the  plate 
circuit  of  the  "normal ly-off"  tube*   The  relay  contacts 
send  plus  or  minus  100  volts  through  the  aileron,  rudder, 
or  elevator  function  switch,  whichever  is  closedo 
The  time  that  the  relay  remains  closed  is  governed  by 
the  RC  time  constant  controlled  by  the  variable  potentio- 
raetero   A  range  of  from  approximately  Oo5  to  2o$ 
seconds  was  availableo   The  electrical  schematic  is 
presented  as  Figo  ?»   A  photograph  of  the  "bread-board" 
model  appears  as  Figo  6„      The  contacts  of  the  relay 
carrying  plus  or  minus  100  volts  would  appear  electri- 
cally just  ahead  of  the  function  switches  in  the  REAC 
solution  schematic  which  appears  as  Figo  9o 

The  six  equations  of  motion  were  put  in  computer 
form  in  Appendix  IIo   These  equations  are  11-13  through 
II-lSo   A  computer  schematic  for  the  simultaneous 
solution  of  these  equation  was  laid-out  and  appears 
as  Figo  9o   From  the  schematic,  patchboards  were  wired 
for  the  lateral  motion  and  control  computer  and  the 
longitudinal  motion  computer »   This  division  of  the 
problem  was  based  on  the  wiring  scheme  v»riiich  would  use 
the  fewest  amplifiers »   The  apparent  symmetry  was 
incidentalo   The  numerical  values  of  the  quantities 
entering  the  problem  were  translated  into  potentiometer 
settingSo  The  values  recorded  were  u  ,  109,   ur  , 
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and  io<^on  one  recorder;  -y  $    lOc^jiOV^,  and  9  on  the 
other  recorder„   The  plus  and  minus  100  volts  was 
placed  on  one  side  of  the  timer  relay  contacts.   The 
other  side  of  the  relay  contacts  was  connected  to  each 
of  the  control  function  switches,  plus  or  minus,  which- 
ever was  required  by  the  schematic* 

The  graph  representing  the  estimated  mass  dis- 
tribution changes  as  the  fuselage  becomes  longer  and 
more  dense  is  presented  in  Figo  lOo   Values  of  the 
moments  of  inertia  of  three  configurations  were  picked 
off  this  grapho   New  stability  derivatives  were  cal- 
culated as  necessary  for  these  configurations.   All  of 
the  moment  derivatives  such  as  Lv  were  affected.   These 
calculated  values  were  translated  into  potentiometer 
settings  and  appear  as  Tables  VII  and  VIII, 

The  configurations  investigated  were  called 
Case  I,  II,  III  and  IVo   The  lateral  computer  was 
referred  to  as  the  **blue"  computer  and  the  longitudinal 
computer  as  the  "red"  computer.  These  are  occasionally 
referred  to  in  the  tables  using  the  abbreviations  B  and 
R. 

The  computers  were  balanced  and  the  recorders 
calibrated.   The  patchboards  were  then  installed  on 
the  computers,   A  rudder,  elevator,  or  aileron  control 
movement  was  then  selected  and  the  recorder  tape  set  in 
motion  at  the  desired  speed.   Pushing  the  time  switch 
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button  sent  an  accurately  timed  standard  step-function 
through  the  control  switch  selectedo   The  solutions  of 
the  equations  of  motion  in  response  to  a  step-function 
input  were  recorded  as  time  histories  of  the  variables 
on  eight  recorder  channels  simultaneouslyo   The  next 
case  was  then  entered  on  the  potentiometers  and  the 
solutions  obtained  as  before*   These  solutions  are 
presented  as  FigSo  12A,  12B,  13A,  13B,  I4A,  and  I4B0 

The  dimensional  stability  derivatives  were  recal- 
culated using  data  for  40,000  feet  and  the  required  angle 
of  attack  increase  obtained  from  Table  IIo  Where 
specific  non-dimensional  derivatives  were  not  available 
for  this  altitude,  those  for  low  altitude  were  extended 
by  using  the  values  of  /^  and  IJ"  at  40,000  feet  in  the 
formulae  of  Tables  III  and  IVo  These  stability  deriva- 
tives and  other  factors  of  the  problem  Wiich  changed  at 
40,000  feet  were  translated  into  potentiometer  settings 
and  run  off  as  before*  The  settings  appear  in  Tables  IX 
and  Xo   The  results  for  Cases  I  and  IV  are  presented  in 
FigSo  15A,  15B,  16A,  I6B,  17A  and  17Bo   It  should  be 
pointed  out  that  the  airplane  which  was  marginal  in 
dynamic  directional  stability  at  sea-level  was  unstable 
at  40,000  feeto   In  order  to  "fly"  even  marginally  at 
this  altitude  the  N^ potentiometer  was  set  to  O0OII5 
and  the  L^  potentiometer  to  zero«   This  had  the  effect 
of  increasing  the  "weather-ccck"  stability  parameter 
about  six  times  and  reducing  the  effective  dihedral 
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parameter  to  zeroo  The  computer  procedure  was  as 
before o   It  should  be  noted  that  solutions  were 
obtained  only  for  Cases  I  and  IV  modified  as  des- 
cribed above o 

The  effect  of  inclination  of  the  longitudinal 
principal  axis  to  the  flight  path  was  investigated 
by  calculating  values  of  moments  and  products  of 
inertis  for  inclination  angles  ^    of  plus  ten,  zero, 
and  minus  ten  degrees o   The  equations  below  obtained 
from  Refo  6  were  useds 

It  was  found  that  only  the  i^^  value  varied  signif icantlyo 
New  values  of  the  coefficients  containing  X^^  were 
calculated  and  translated  into  potentiometer  settingSo 
These  are  presented  in  Table  XIo   Several  sign  changes 
were  requiredo   These  were  accomplished  by  using  inverters 
as  necessaryo   The  computer  was  operated  as  before 
using  Case  I  at  low  altitude* 

The  C^  and  C.     programs  were  accomplished  by 
changing  these  values  over  the  range  considered, 
recalculating  the  affected  stability  derivatives,  and 
translating  them  into  potentiometer  settingSo   This 
was  done  for  Case  I  at  sea  level  and  at  40,000  feeto 
These  values  are  presented  in  Tables  XII  and  XIIIo 
The  computer  operation  was  similar  to  the  other  pro grams o 
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RESULTS  AND  DISCUSSION 

Time  histories  of  the  variables  U  ,  ur-  ?  0  >  tt  , 

* 
<p      i    Q     9   f  9    and  >//  were  recorded  as  the  airplane 

responded  to  a  step  function  control  input  for  each 
change  in  a  parameter  being  investigatedo   To  study 
the  trend  of  the  effect  of  the  parameter  changes,  the 
curves  were  superimposed,  each  curve  being  drawn  with 
a  characteristic  lineo   These  composite  solution  curves 
are  presented  as  FigSo  12A  through  23Bo   The  standard 
step  function  time  plot  appears  as  Figo  llo   It  should 
be  noted  that  the  A  and  B  parts  of  the  figures  each 
represent  one  of  the  four  channels  which  were  simultan- 
eously recordedo  Where  separate  lines  for  each  para- 
meter change  can  not  be  observed,  it  should  be  assumed 
that  the  curve  lies  behind  the  previous  curve,  or  curves, 
and  that  no  significant  change  was  caused  by  the  new 
value  of  the  parameter o 

The  mass  distribution  program  is  presented  in  the 
FigSo  12A,  12B,  13A,  13B,  14A,  and  I4B0   The  elevator 
response  curves  show  no  significant  changes  caused  by 
the  changed  moments  of  inertia o   The  elevator  set-up 
the  phugoid  motion  vdiich  reriHined  essentially  constanto 
The  short  period  motion  was  virtually  dead-beat  dampedo 
The  amplitude  of  Q   in  Case  IV  is  reduced  since  the  air== 
craft  was  more  sluggish  to  pitch  change  due  to  the  in- 
creased pitching  moment  of  inertiao   A  slight  dutch- 
roll  was  inducedo   This  was  not  anticipated  since  the 
equations  indicated  that  there  would  be  no  coupling 


from  the  symmetric  motion  into  the  lateral  motion  if 
there  was  no  initial  lateral  motion,,   A  possible 
explanation  might  be  that  there  was  some  unbalance 
in  the  servo-multipliers  or  jitter  in  the  lateral 
motion  computer  due  to  small  line  current  variationso 
The  same  effect  might  be  observed  in  flight  tests 
caused  by  slightly  unsymmetrical  control  application 
or  small  trim  unbalanceo   Spiral  instability  was  set- 
up presumably  by  the  same  mechanism  which  caused  the 
slight  dutch-roll  to  start » 

The  cross-coupling  set-up  by  the  aileron  and 
rudder  motion  can  be  observed  in  the  or  ploto   The 
coupling  caused  by  the  more  assymmetric  motion  induced 
by  the  rudder  step  is  several  times  that  of  the  aileron 
cross-couplingo  The  aileron^s  cross-coupling  decreased 
as  the  values  of  i^^  and  T-^  were  increasedo   This  was 
unexpected  but  can  possibly  be  explained  by  the  fact 
that  two  opposing  effects  were  taking  placeo  The 
coupling  increased  but  the  airplane  became  more  slug- 
gish since  the  pitch  and  yaw  moments  of  inertia  were 
increasedo   This  sluggishness  can  be  observed  in  the 
progressive  decrease  in  amplitude  and  frequency  of 
and  T  o   The  large  decrease  in  roll  moment  of  inertia 
had  very  little  effect  on  '^  and  <^  o 

The  altitude  effects  on  coupling  and  dynamic 
stability  are  presented  in  FigSo  I5A5  15B,  I6A5  I6B, 
17Aj  and  I7B0  Many  of  the  observations  made  in  the 
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sea-level  program  apply  here.   The  short  period 
oscillation  damping  decreased  increasing  the  period 
of  this  motion*   The  aileron  coupling  was  larger  for 
Case  IV,  that  is,  as  the  fuselage  became  long  and  dense. 
The  same  effect  was  observed  in  the  coupling  caused  by 
the  ruddero   This  time  there  was  little  amplitude 
difference  between  the  aileron  and  rudder  coupling. 
The  sluggishness  caused  by  increased  I^  and  I^  , 
together  with  decreased  aerodynamic  damping,  made  the 
swings  progressively  larger  in  amplitude  and  period. 
The  roll  angle  did  not  significantly  change. 

In  the  program  which  considered  the  effect  of 
changing  the  inclination  of  the  longitudinal  principal 
axis  to  the  flight  path,  the  Ij^^  program  essentially, 
only  the  rudder  and  aileron  step  input  response  was 
investigatedo   As  Yj  became  more  positive,  dynamic 
stability  increased;  more  negative,  dynamic  stability 
decreased.   At  an  Vj  of  minus  ten  degrees,  dynamic  in- 
stability  developed.   From  the  coupling  standpoint,  the 
same  effect  was  observed.   There  is  less  coupling  when 
V(   becomes  more  positive;  more  when  ^  becomes  more 
negative.  The  rudder  curves  show  these  effects  more 
clearly  than  the  aileron  curves. 

The  effects  caused  by  varying  the  "weather- cock" 
parameter,  Cp,  ,  is  presented  in  Figs,  20A,  20B,  21A 
and  21Bo   Only  the  aileron  and  rudder  farcing  functions 
were  considered.  The  curves  for  two  values  of    ^n» 
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at  high  altitude  were  superimposed  on  the  sea-level 
curveso   At  both  altitudes  similar  effects  were 
observedo   Coupling  increased  as  C   decreasedo   The 
dynamic  lateral  stability  decreased  as   C   decreasedo 
This  was  to  be  expected,  of  course,  since  decreased 
tail  arm,  or  tail  area,  lead  toward  lateral  instabilityo 
As  the  angular  velocity  changes  increased  with  this 
instability  the  coupling  increasedo 

The  curves  for  the  Cj^  program  are  presented 
in  FigSo  22A,  22B,  23A  and  23Bo   The  entire  sea^level 
program  was  presented  but  only  the  curve  for  a   C^ 
of  0<,20  at  40,000  feet  was  addedo   As  C.  was  increased 
dynamic,  lateral  instability  quickly  developed  along 
with  increasing  couplingo  The  effect  was  more  apparent 
in  the  motion  caused  by  the  rudder  forcing  functiono 
The  reduced  damping  at  altitude  hastened  the  develop- 
ment of  instability  and  the  coupling  increasedo   The 
curve  for  a  Co     of  Oo20  was  marginally  stable  at  sea 
level  but  unstable  at  40,000  feeto 

The  parameters  C^  and  Cj^  had  very  significant 
effects  on  dynamic  stability*  The  effective  dihedral 
parameter,  Cj^  ,  was  very  criticalo  This  would  appear 
to  make  analysis  of  the  dynamics  of  the  airframe  very 
important  in  the  design  stage  of  swept-wing  and  delta- 
wing  air  craft o  The  REAC  provides  a  powerful  tool  for 
such  an  analysis o 
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The  use   of   a  standard   step-function   input 
provided   a   good  means  for   studying  the   effect   of  a 
parameter   changeo      However,    this  method  was   inadequate 
from  the    standpoint   of   bringing   out   any   instabilities 
caused  by  cross-coupling  of  the   type  investigated   in 
Refo    4o 

It  should  be   pointed  out  that  the    dynamic  be- 
havior observed  was  that   of  one   airplane ,    or  variation 
of  one   airplaneo      Although  effort  was  made  to  make 
the  mathematical  model   investigated   "typical"    ,    it 
would   not   be   reasonable  to  suggest   that  all  airplanes 
behave   in  the   manner   of  the    airplane  investigatedo 
On  the   other  hand,    the   procedures  used   in  the   analysis 
were  straight-forward  and   it  is    felt  that  the  dynamic 
behavior  observed  might   be   considered  reasonably 
typicalo 
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CONCLUSIONS  AND  RECOMMENDATIONS 

1.   Mass  distribution  changes  had  no  significant 
effect  on  the  motion  of  the  aircraft  in  the 
longitudinal  plane. 

2o   The  motion  of  the  aircraft  in  the  symmetrical 
plane  coupled  into  the  lateral  motion  in  the 
form  of  a  slight  dutch  roll  and  also  triggered 
the  spiral  instability  raodCo   However,  this  is 
believed  to  have  been  due  to  slight  assymmetry 
in  the  computing  systemo 

3<»   The  assymmetrical  motion  of  the  airplane  coupled 
into  the  motion  in  the  symmetrical  plane  and  is 
most  significant  as  a  change  in  u/  »   This  would 
be  equivalent  to  a  change  in  the  angle  of  attacko 

4«  The  coupling  effect  increased  as  the  fuselage 
became  longer  and  more  dense. 

5.  The  assymmetrical  motion  caused  by  the  rudder 
had  a  greater  coupling  effect  than  the  motion 
caused  by  the  ailerono 

6.  At  high  altitudes  both  the  coupling  and  the 
trend  toward  instability  increased. 

7»   At  a  positive  inclination  of  the  principal 
longitudinal  axis  to  the  flight  path,  there 
was  less  coupling  and  increased  dynamic  stability* 
A  negative  inclination  angle  had  the  opposite 
effect o 
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So  A  decrease  in  the  "weather-cock"  stability 
parameter,  C^  ^  resulted  in  a  decrease  in 
dynamic  stability  and  an  increase  in  couplingo 

9.   An  increase  in  the  effective  dihedral  parameter, 
Cji^    ,  resulted  in  a  decrease  in  dynamic 
stability  and  an  increase  in  couplingo 
lOo   The  standard  step-function  input  which  supplied 

the  forcing  function  provided  a  means  for  studying 
the  effect  of  a  parameter  change  on  the  coupling 
and  the  dynamic  stability  of  the  airplaneo   How- 
ever, it  does  not  appear  adaptable  to  the  invest!  ■= 
gat ion  of  instabilities  which  may  develop  from 
the  coupling  caused  by  rapid  rolling  or  yawingo 
An  investigation  might  do  this  by  using  a  pro- 
grammed input  which  would  reverse  before  the  REAC 
limits  were  reachedo 

The  dynamic  behavior  observed  in  this  investigation 
was  that  of  the  mathematical  model  of  one  airplane  or 
variation  of  one  airplaneo   Effort  was  made  to  use  a 
mathematical  model  which  could  be  considered  reason- 
ably "typical "o   However,  it  would  not  be  reasonable 
to  suggest  that  all  airplanes  wDuld  respond  similarlyo 
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SYlvlBOLS 

Uo  initial  airspeed,  feet  per  second 

0  mass  density  of  air,  slugs  per  cubic  foot 

b  vd.ng  span,  feet 

S  wind  area,  square  feet 

W  weight  of  airplane,  pounds 

rf\  mass ,  slugs 

g  acceleration  of  gravity,  feet  per  second  per  second 

AR  aspect  ratio 

M  Mach  number 

C^  wing  chord  tip,  feet 

C^  wing  chord  root,  feet 

I^  moment  of  inertia  about  principal  longitudinal 
axis,  slug-feet^ 

Xy  moment  of  inertia  about  principal  tranverse 

2 
axis,  slug-feet 

I   moment  of  inertia  about  principal  normal  axis, 

slug-feet^ 
Zj^^  moment  of  inertia  about  flight  path  axis,  slug-feet 
I^  moment  of  inertia  about  the  pitching  axis, 

slug-f eet^ 

Ij    moment   of  inertia  about  the  yav^^ing  axis,    slug-feet^ 
^xi  P^<^duct   of   inertia  with  respect  to  the   flight   path 
axis  and  axis  normal  to  flight  path,    slug-feet^ 
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p   force,  pounds 


h       moment  of  momentum,  (slug-ft)-ft  per  second 

C   lift  coefficient  ,  -  '"'^^  ^"^^ 

rolling  moment  coefficient  ,  2 _ 

i  ^  '      Cj  sb 

C,  lateral-force  coefficient  ,     lateral  fer-ce 

Q^    aerodynamic  pitching  moment  coefficient  about  the  e.g. 

Cn  drag  coefficient  ,   Prae  fo.-cc 

<^   angle  of  bank,  radians 

Q      angle  of  pitch,  radians 

1^   angle  of  yaw,  radians 

S       control  deflection,  radians 

S       angle  of  sideslip,  radians 

y   flight  path  angle;  angle  between  relative  wind 
and  the  horizontal,  radians 

y^       angle  between  the  longitudinal  axis  of  the  fuse- 
lage and  the  relative  wind,  radians 

F   angle  between  the  thrust  line  and  the  longitudinal 
axis  of  the  fuselage,  radians 

2'   perpendicular  distance  from  eg,  to  the  thrust 
line,  feet 

y   fraction  of  critical  damping 

CO   frequency,  radians/sec. 

\j^      initial  forward  velocity  in  direction  of  positive 
axis  (feet  per  second) 

■\/"   initial  side  velocity  (feet  per  second) 

W"  initial  normal  velocity  (feet  per  second) 
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U  forward  velocity  change  after  disturbance  (feet 

per  second) 
u  side  velocity  change  after  disturbance  (feet  per 

second) 
ur  normal  velocity  change  after  disturbance  (feet 

per  second) 
P  initial  rolling  angular  velocity  (radians  per 

second) 
Q  initial  pitching  angular  velocity  (radians  per 

second) 
R  initial  yawing  angular  velocity  (radians  per 

second) 
•p^^  rolling  angular  velocity  increment,  radians  per 

second 
a^e   pitching  angular  velocity  increment,  radians 

per  second 
/L^^/  yawing  angular  velocity  increment,  radians  per 

second 
Cr)  the  change  in  drag  coefficient  with  varying 

u. 

forward  velocity 
Q       the  change   in  drag  coefficient  with    varying  angle 

of  attack,    per  radian 
C^     rate   of  change  of  lateral-force  coefficient  with 

rolling-angular-velocity  factor,    per  radian 
O     rate  of  change   of  lateral-force   coefficient  with 

yawing-angular-velocity  factor  ,    per  radian 
C^     lateral  force  derivative,    rate    of  change   of  lateral 

force   coefficient  with  angle   of  sideslip,    per 

radian 
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Cy   the  change  in  side  force  coefficient  vdth  variation 

in  rudder  deflection,  per  radian 
Cy       the  change  in  side  force  coefficient  with  aileron 

deflection,  per  radian 
C^      change  in  lift  coefficient  with  variation  in  forward 

velocity,  angle  of  attack  and  attitude  remaining 

constant 
C^   change  in  lift  coefficient  with  varying  pitching 

velocity,  no  change  of  airplane  angle  of  attack 

as  a  whole ,  per  radian 
C,   the  change  in  lift  coefficient  vdth  varying  angle 

of  attack,  "lift  curve  slope",  per  radian 
C^   change  in  lift  coefficient  with  changes  in  elevator 

deflection,  per  radian 
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Q       rate  of  change  of  rolling  moment  coefficient  with 
yawing  angular  velocity  factor,  per  radian 

Q  effective-dihedral  derivative,  rate  of  change  of 
rolling  moment  coefficient  with  angle  of  side- 
slip, per  radian 

C.      damping  in  roll  derivative,  rate  of  change  of  rolling 
moment  coefficient  with  rolling-angular- velocity 
factor,  per  radian 

C   change  in  rolling  moment  coefficient  with  variation 
in  rudder  deflection,  per  radian 

Q   aileron  effectiveness  on  aileron  power,  the  change 
in  rolling  moment  coefficient  with  change  in 
aileron  deflection,  per  radian 
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C^   longitudinal  static  stability  derivative,  change 

in  pitching  moment  coefficient  with  varying 

angle  of  attack,  per  radian 
C^.  the  change  in  pitching  moment  coefficient  l^d.th 

variation  in  rate  of  change  of  angle  of  attack, 

per  radian 
C^   pitch  damping  derivative,  change  in  pitching  moment 

coefficient  with  varying  pitch  velocity,  per  radian 


"^^ 


Crr^      elevator  effectiveness  or   elevator  power,   the  change 
in  pitching  moment   coefficient  with    change  in 
elevator  deflection,    per  radian 

C„      directional   stability  derivative,    weather-cock 

Stability,    rate  of  change   of  yawing-moraent   co-= 
efficient  with  angle   of   sideslip     \ -r-^ )     per  radian 

C„      damping  in  yaw  derivative,    rate   of   change   of   yawing 
moment  with  ya\'d.ng  angular  velocity  factor. 


per  radian     (  TTb~l 

■5 1  r.  / 


•  n 


On   rate  of  change  of  yawing  moment  coefficient  with 

rolling  angular  velocity  factor,  per  radian  (  ^^ -p b 
C    rudder  effectiveness  or  rudder  power,  the  change 

in  yawing  moment  coefficient  with  variation  in 

rudder  deflection,  per  radian 
C^   change  in  yawing  moment  coefficient  with  change  in 

aileron  deflection,  per  radian 
Xu   change  in  fore  and  aft  force  with  change  in  forward 

velocity,  per  second 

Xq   change  in  fore  and  aft  force,  change  in  pitching 
I* 

velocity,  feet  per  second-radian 
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y,^     change   in  fore  and  aft   force   vdth   change  in 

normal  velocity,    per   second 
Xj^    change   in  fore  and  aft  force  with  rate  of   change 

in  normal  velocity,    per   foot 
X^     change   in   side   force   with  change   in  yawing  velocity, 

per   radia  n 
Y^      change   in  side   force  mth   change   in   side  velocity, 

per   second-radian 
Y-p      change  in   side  force  with  change   in  rolling  velocity, 

per   radian 
Yg      change  in  side  force  with    aileron  deflection, 

per  second=-radian 
V^      change   in  side  force  with   rudder  deflection,    per 

second-radian 
Zl^     change   in  normal   force  with  varying  forward  velocity, 

per   second 
Z<i     change   in  normal  force  with  varying  pitching  velocity 

feet  per  second-radian 
"Z       change  in  normal  force  with  varying  vertical 

velocity,    per  second 
E^^     change   in  normal  force  with  elevator  deflection, 

feet   per  seconds-radian 
L       change  in  rolling  moment  vdth   varying  yawing  velocity 

per   second-radian 
L        change  in  rolling  moment  v/ith  varying  side  velocity, 

per  second-feet 
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L  p  change   in  rolling  moment  with  varying  rolling 

velocity,    per   second-radian 
L^  change   in  rolling  moment  with  rudder  deflection, 

per  seconds-radian 
Lg    change   in  rolling  moment  with  aileron  deflection, 

per  seconds-radian 
My  change   in  pitching  moment   with   varying  forward 

velocity,    per  second-feet 
M^r  change   in  pitching  moment  with  varying  vertical 

velocity,    per   second-feet 
M^  change   in  pitching  moment  with   rate   of  change   of 

vertical  velocity,    per  foot 
Mq    change   in  pitching  moment  with  varying  pitching 

velocity,    per   second-radian 

M^  change  in  pitching  velocity  with  elevator  deflection, 

2 
per   second  -radian 

My.  change   in  yawing  moment  with  varying  side  velocity, 

per  second-feet 
Klj^    change   in  yawing  moment  with  varying  yawing 

velocity,    per  second-radian 
M      change   in  yawing  moment   with  varying  rolling 

velocity,    per  second-radian 
Ng    change   in  yawing  moment  with   rudder  deflection, 

per  seconds-radian 
Ng    change   in  yawing  moment  with   aileron  deflection, 

per  s ec ond2-=  radian 
M     change   in  yawing  moment  with  side-slip  angle, 

per  second  -radian 
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(  )o  refers  to  initial  steady  flight  condition 

(  )^  pertains  to  the  fore  and  aft  axis 

(^  )  pertains  to  the  transverse  axis 

(  )^  pertains  to  the  normal  axis 

C  )f^  elevator 

(^  ]g  aileron 

C  ).  rudder 
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APPENDIX  I 

Development  of  the  Equations  of  Motion 

The  equations  of  motion  of  the  airframe  are  written 
by  equating  the  forces  and  moments  acting  on  the  air- 
frame to  the  planers  reactions  in  accordance  with 
Ne\>rton^s  lawso   The  reference  axes  are  considered 
fixed  in  space. 

The  aircraft  is  considered  to  be  a  rigid  body; 
aeroelastic  effects  are  not  considered.   The  motion 
is  referred  to  right-handed  Cartesian  axes.   The  axes 
used  and  the  sign  conventions  observed  are  shown  in 
Fig,  3  and  Table  lo 

Newton's  second  law  states  that  the  rate  of  change 
of  momentum  of  a  body  is  proportional  to  the  force 
applied  to  the  body  and  the  rate  of  change  of  the 
moment  of  momentum  is  proportional  to  the  torque  applied 
to  the  bodyo   The  mass  of  the  airplane  is  considered 
constant  for  the  duration  of  the  dynamic  maneuver 
considered, 

U)  I  L  -  dh^ 

at 

dt 
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Investigating  the  dynamics  of  an  infinitesimal 
element  of  mass  dm      of  thie  airplane  referred  to  the 
Cartesian  coordinate  system  yields  equation  set  (3)o 

rfhy    •     (iSx')Qdrr,    -   x-^    Pdm     -L^^     R.clr/1 

For  a  finite  mass  these  become 

(4.)  U^  =     pfcci2  +  2^)dm    -    Q|x<^dm    -  Rfxz-dr 

The    integral         (n^-t-a*-)  dm   is  recognized  as  the 
moment   of   inertia     I^    of   the  whole  airplane  about   the 
X-axis o      Equation  set    (4)    can  be  written 


Im 
irn 


K- 

PI, 

-  O.  Uy 

-K  r,^ 

s- 

Oly 

-  "2  II* 

-    Pl.y 

The  equations  of  motion  relative  to  inertial 
axes  becomes: 
^  5  J  ^  dt 

I  f  y   =     '^    lY 

dt 

d-t 

Zm^  ^^  6iy  t  Qi^  -  R  I,,,  -  Ri\^-  Pix^ '  Pi^x^t 

Equations  (6)  are  written  with  respect  to  fixed 
inertial  axesc   It  is  expedient  to  use  Eulerian  axes 
with  their  origin  at  the  c.3.  which  are  fixed  in  the 
airplane  and  move  with  ito   At  any  instant  the  airplane 
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in  a  maneuver  has  a  motion  referred  to  the  Eulerian 
axes  and  the  axes  have  a  motion  referred  to  fixed 
spaceo   The  earth  and  its  atmosphere  are  regarded  as 
fixed  space.   The  absolute  acceleration  of  the  aircraft 
referred  to  the  Eulerian  axes  can  be  written 

where  V^  is  the  total  velocityo 

Similarly  the  absolute  rate  of  change  of  the 
moment  of  momentum  is  written 

(B)  Ih^^*   d_b  +  z:*^  h 

dt      dt 
The  first  term  arises  from  the  change  relative  to 

the  Eulerian  axes  fixed  from  instant  to  instant  in  the 

airplane  and  the  second  considers  the  change  of  the 

Eulerian  axes  relative  to  inertial  axes  fixed  in  spaceo 

Moments  and  products  of  inertia  measured  relative 
to  Eulerian  axes  are  independent  of  time  since  the  mass 
is  assumed  constant  for  the  duration  of  a  maneuvero 
Many  of  the  terms  in  equation  set  (6)  go  to  zero  for 
this  reasono 

The  term  loxV/^  of  equation  (7)  can  be  written: 

(Sa)  "tiycs/^  ^      i     j   L 

P  Q  R 
U  V  ^  . 

Where  il  ^  j  ^  k  are  unit  vectors  directed  along  the 
Eulerian  X,  Y^  2.  axes. 

Expanding  the  cross  product  and  combining  with  the 
components  of  the  first  term  of  (7)  that  equation  can 
be  written  in  component  form  as: 


d,  r  ^  +  pv-  au 
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By  parallel  reasoning    (8  )    can   be  written^ 
(10)     ^'L/    ^^  i^Q-S*^^  PL -o^-il,, -QRU,-I,) 

-POTx*  -q"I?.|  -M^Mya     -^PJ^I,^ 


-QR  T<^  -  1^*  I^a.  +  P*"  lyfr    -»-   <^a  I 


*i 


dh, 


With  the  Eulerian  axes  as   chosen  and  shown  in 
Fig.    ^,    the     Kz   -plane   is  a   plane   of    symmetry.      There 
is  both  a    positive  and   a   negative  value   of    y     for 
each  value  of   z     ;    consequently     Ty^  ■  (tjtdw^o  and 
1^^  -     J  >^^  ^^  =   O  . 

The   equations   of  motion  of  an  aircraft  referred  to 
Eulerian  axes   can   then  be  written: 
(11)         ?F^  -  /VV1  [l/ +  Qv^/-RVJ 

Z  Ft    -  '^nC'A/  +P\/  -QUj 

IL-     PI,  -  R  Ix-t   -^  QRCli-Iv)  -  PQl^t 

InJ     -     Rl^  '  Pl/t  +  PQ(I<^-I>)  ^  Q^Iyt 
The   left   sides   of  the  equations   in  set    (11)   are 
the  summations   of  the  external  forces  and  moments 
applied  to  the  airplane   in  flight.      These   external 
forces  are  gravity  forces,    aerodynamic  forces,    and 
thrust  forces. 

The  gravity  force  may  be  considered  to  act  at  the 
center  of   gravity  making  no   contribution  to  the 
summation  of  moment So 
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The  resolution  of  the  gravity  force  into  components 
referred  to  disturbed  Eulerian  axes  is  shovm  in  Figo  lo 

The  force  relations  may  be  written 
(12)  IF^'=  rvvi[U+QW  -  RV]  -  X3 

iFj'  -   /^[vv  ^PV-  QU]  -  l^ 
by  transferring  the  gravity  terras  to  the  right-hand 
side  of  the  equationo   The  gravity  terms  /^ ,  Y3  ,  and 
Z^     are  expanded  in  Figo  1» 

The  left  side  of  all  the  force  and  moment  equations 
now  include  aerodynamic  and  thrust  forces  as  well  as 
moments  due  to  control  surface  deflection. 

The  disturbed  motion  of  an  airplane  at  any  instant 
can  be  considered  the  result  of  disturbing  the  airplane 
from  some  steady  flight  condition.  This  motion  may 
be  referred  to  as  disturbed  flight  and  the  Eulerian 
axes  as  the  disturbed  axes.  The  forces  acting  on  the 
airplane  .under  these  conditions  may  be  regarded  as 
the  sum  of  the  force  required  to  produce  steady  flight 
plus  an  increment  required  for  the  disturbance.  Accord- 
ingly, each  of  the  total  instantaneous  velocity  com- 
ponents can  be  written  as  the  sum  of  a  constant  velocity 
component  during  steady  flight  and  a  change  caused  by 
the  disturbance: 
(13)  U  =  Uo  -^^ 


Substituting  in  (12)  considering  that  ^  =  i3. ,  o 

dt      dt         ' 

etCo,    and   using  the    expansion  of  the   gravity  terms 
from  Figo    1  after   setting   terms   like     s\n4^s[r^4•^  o     , 
5irvv|y-v4/      and      cosG--!     ,    the  set   of   equations   below 
can  be  written o      The   assumptions  that  the    sine   of  an 
angle  can  be  replaced   by  the   angle,    the   cosine   is 
equal  to   one  and    the    product   of  sines   can  be  replaced 
by   zero   produces  results  vhich,    although  approximate, 
are   in  good  agreement  v/ith  physical  results   according 
to  Refo    7o      These   assumptions,    of   course,   apply  to  the 
disturbed  angleSo 

^(Po^t      ZPolP    +-^0     lyt 

IM  ^  AU  -t  Ixi  ^  (f'oQc  -^^o^  +Qo>  r>^l(l,-rx)  ^(OoR,-^Qo^  ^(?o^-^a(^)l^^ 

The  aerodynamic  forces  and  moments  acting  on  an 
airplane,    the   left   side   of  the    above   equations,    are 
functions   of  the   flight    conditions   and    of  the   deflection 
of  the   control  surfaceso 
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Each  of  the  forces  {/,N',Qnd2.  )  and  the  moments 
(  L_M,<J'^ci^)  )  can  be  expressed  for  small  oscillations 
by  expanding  by  MacLaurin^s  theorem  neglecting  second 
order  terras  to  obtain: 

where   IX  ,  1>^   etc.  are  partial  derivatives  of  X 
with  respect  to  w,q^j  etc.   Exactly  analogous 
expansions  are  made  for  V,  2,  L,  f^i,  and  f^  o  These 
forces  and  moments  can  be  expressed: 
(15)  K-  Xl-  H^-  l^a-^M-^  *  M.  ^  •  -  '  • 


^'^     au-      d^^  f^    ^v     ^Sn.  '^ 

Each  of  the  terms  in  (15)  has  a  physical  sig- 
nificance,  X'o  Yl  and  2„   L  Mo  and  ^c    are  the  forces 
and  moments  acting  along  and  about  the   x  >  H,  Q'^'^'  Z: 
axes  \^ile  the  plane  is  in  steady  flight*   The  terms 
similar  to  ^  u.  express  the  change  in  the  given 
force  or  moment  caused  by  the  given  disturbance  quan- 
tityo   The  coefficients  of  the  form   ^u   are  the 
basic  dimensional  stability  derivatives.,  A  complete 
discussion  of  stability  derivatives  may  be  found  in 
Refs.  (1),  (2),  (3),  (7),  (3),  and  (9). 


i 
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Thrust  effects  are  part   of   the   left   sides  of 
equation  set    (14)o     The  equations   for  steady  flight 
condition  become: 
(16)  y.  =  t;  CO.  I 

Mo-     X^j 
Since  the  Eulerian  axes  remain  fixed  vri.th   reference 
to  the   airplane  during  a  disturbance,    the  thrust 
components  relative  to  disturbed  axes  become: 

v^ere    T|   the  thrust  during  the   disturbance,    is    T^  •♦-  aT« 
AT~  can  be   considered 

(18)  AT-    ^-   ^^    S,,^ 

and 


i  X         fH'--^ 


(19)  X-T,  cost    ^  tc^'S)lT^    ^(^cc^pjiTS 

•Jut  ^    ^^ni^^ 

The  equations  for  steady  flight  are  found  by 
substituting  the  expressions  obtained  above  for  the 
thrust  and  the  aerodynamic  forces  (and  moments)  into 
the  left  sides  of  equation  set  (14)  and  setting  the 
disturbed  quantities  equal  to  zero: 

(20)  X,  -  w^m  ©^  -t-  "n,  Cos^  -  O 

Yo-^  o  i-  o       -a 


-^4- 

Lo  +  o  -t-  o  *   o 

Mo  -^  O  +  ITij  -  O 

kl^  r    o  t    o     -  o 

If  now  the  complete  equations  for  the  disturbed 
airplane  are  written  with  both  steady  flight  and 
disturbed  terms  they  will  be  of  the  forms 
Aerodynamic  forces  +  thrust  forces  =  inertia  force  <= 
gravity  force 

Equation  set  (20)  can  then  be  subtracted  from  the 
complete  equations  of  motion  for  the  disturbed  air- 
plane yielding: 

-  P0VV9  -Pour  -  l^/o-j?-ur>  ~  Cssu^^ol'f  -300^0  sw,  (|)^  -  (gco^Go  c<r3^J<^J 
=  ^»^  [or  +  Pa  V«   +-  PoT  +v/oi)  -t-T'>r    -  OoTJ;  -QoU^-TTo^^-C^cL 

-^  I.^  +CPoRo+H)L  ■^-i?o>^■>rL)(Iy-It) 


y* 


■35- 


Dividing  the  force  equations  by  rwi  and  the  moment 
equations  by  the  appropriate  moment  of  inertia  yields 
terms  of  the  form: 

Replacing  —  Z>y  »  by  Vu.  and  J_  :)L  n/  by   L,   simplifies 
the  notation.   These  quantities  are  called  dimensional 
stability  derivatives,   A  list  of  these  stability 
derivatives  in  terms  of  the  dimensionless  forms  appears 
as  Tables  III  and  IV.   The  flight  path  angle  l^  replaces 
Oo  •    The  equations  in  final  form  can  be  written; 
(22}  %t^^-^^^%  ^"/u.ur^^dr^   ^G^J|  T^^L  4-Qos^  Tj^  S^ip^ 

=   U.   -t-     W„<^  •^  Oour  -\  uJ-e^'  Q^\/^  -|2.(5Tr    -  ^4  JL    -  '\rK. 
■•-'^a  +    Urt.  -  Pour  -(c,  3u^>lo)f  -  g  5.:^.  V,,  Si:^  4^,  -  C j  CioSo    c«34'o)  4 

»     kIt^T-H  >viR./L4  a^-pH  K)g    S^    -^KJ        S, 


-JL    ■■  "       "^ 


APPENDU  II 
Adaptation  of  the  Equations  of  Motion 
to  the  Analogue  Computer 

The  dynamic  equations  of  motion  of  the  airplane 
were  developed  in  Appendix  lo   These  are  summarized 
below: 

(ir-s)  ^w.^-+  2<^^-f  ^^ur  +  "2:s^Sj_~  sv^  ^  T:^w  -  ^S^J^,^^v.^  * 

(11-4-)  ^  - ft^ivB-  -v-Cqo'^s  +  QoiL-  f^^  +  ^n.)  C  r^-ii^) 

T  1 

"^  I*- 

These  equations  are  general  and  allow  initial 
velocities,  angular  or  linear,  with  respect  to  any 
axiso   Additional  control  terms ,  such  as  flap 
deflection  may  be  included  by  adding  the  product  of 


i 


-87- 

the  desired  deflection  and  the  appropriate  stability 
derivativeo   Any  attitude  of  the  airplane  is  also  per- 
mittedo 

In  order  to  fit  the  two  analogue  computers  avail- 
able certain  assumptions  were  necessaryo  The  plane  was 
assumed  to  be  in  horizontal  level  flight  with  no  initial 
velocities  except   "IX,  »  forward  velocity  along  the 

)<-  -axiso   Stability  axes  are  chosen  with  the  undis- 
turbed X--axis  alligned  with  the  relative  wind  and  the 
>/  -axis  a  principal  axiSo   The  product  of  inertia  will 
be  different  from  zeroo   The  terms  which  describe  power 
effects  are  permitted  to  go  to  zero  by  assuming  gliding 
flight  since  these  power  terms  are  not  essential  to  the 
dynamic  effects  which  are  studiedo  They  merely  add 
complexityo 

Solving  the  resultant  equations  in  terms  of  the 
highest  derivative  places  them  in  the  form  required  for 
analogue  computer  solution: 
(ir-l)  u  •  Xu.tA  +  ^<^^  -^  ^ur  i^  -+-  >^t>  '^    -t-  "^n-  -  ^a.  -  <g  e 


H.  "       "fl^ 


"^  %    "^t  l^^^  -^^        l^ 


^-  '  I. 


The  assumptions  made  simplify  the  equations 
sufficiently  to  permit  use  of  the  analogue  computerso 
The  inertial  coupling  terms  remaine   These  products 
are  handled  by  the  servo-multipliers  and  permit  the 
airframe  studied  to  move  vd.th  six  degrees  of  freedomo 

The  dimensional  stability  derivatives  are  obtained 
from  Table  V,  substituting  these  values  in  the  equations 
above  yields  the  analogue  computer  equations  for  the 
basic  airframe: 

Longitudinal  Equations 
(jr-  \i)      U  *  — •  oi  5  1  u„  -+  ,  o  -LGH  <^  -+  'VJl.  "  ^^<^  "  3z.Z© 
(^-  14)   Ur  *  -.  130  u.  H-   (I.  ^S"  ^  "  "Z-IO  ur  +  (oLS  Sjj. 

(jr-l5)  4'  -.OOOZSlw.-  IU.649<jj-.0?4MJ--''14..d5j  -.00184  i>  +-.868-^^.  t.05fl^-.05  '^ 

Lateral  Equations 
(jr-ifa)  V- *  +  \.&i  n--.  nsv  -  50t-|3  +Z1  .4-^;^. 

These  equations  of  motion  are  entered  directly 
in  the  analogue  computers «   The  dimensional  form  of 
the  stability  derivatives  permits  the  use  of  real  time, 
that  is,  one  second  of  machine  time  is  equal  to  one 
second  of  real  time.   The  coefficients,  which  are  the 
dimensional  stability  derivatives,  are  entered  as 
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potentiometer  settings  either  directly  or,  if  they 
are  greater  than  one,  as  a  fraction  which  is  sub- 
sequently ajnplified  to  full  valueo   The  potentiometer 
settings  along  with  the  parameters  they  represent 
are  presented  in  Tables  VI  and  VII. 

The  stability  derivatives  which  represent  the 
change  in  a  variable  due  to  a  control  movement  are 
based  on  a  one  radian  control  displacemento   The 
product  of  the  control  movement  and  the  corresponding 
stability  derivative  represents  the  forcing  function 
vdiich  disturbs  the  airframe  from  its  level  flight 
equilibriumo   Since  the  equations  are  entered  directly 
with  no  machine  scaling  there  is  a  one  to  one  corres- 
pondence between  machine  units  of  100  volts  and  unit 
movement,  or  in  the  case  of  angular  control  movement 5 
unit  radian  throwo   For  a  10°  control  movement  a  ratio 
is  set  up  which  establishes  the  potentiometer  setting 
for  a  10^  throwo   A  lesser  or  greater  movement  would, 
of  course,  be  represented  by  a  proportionately  greater 
or  lesser  potentiometer  setting.  An  example  is  calcul- 
ated below  to  obtain  the  potentiometer  setting  represent- 
ing  M^  for  a  10^  control  movement « 
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That  is  to  say  that  v^en  one  machine  unit,  or 
100  volts,  is  placed  across  the  M^  potentiometer  set 
at  .765,  the  analogue  of  the  change  to  the  pitching 
moment  caused  by  an  elevator  throw  of  one  radian  will 
enter  the  equations  of  motion  as  a  forcing  function,, 
Ten  degrees  would  be  this  value  multiplied  by  the 
fraction  -i^   or  .134«   The  other  potentiometer 
settings  representing  the  change  caused  by  control 
movements  are  similarly  computed «   For  the  basic  air"= 
plane  these  values  appear  in  Tables  VI  and  VIIo 

The  computer  diagram  for  the  simultaneous  solution 

of  these  six  non-linear  equations  of  motion  is  presented 

as  Fig,  9o   The  inertial  coupling  terms,  which  are  the 

product  or  moment  of  inertia  factors  and  velocity  terms 

such  as  ^n.    ^J^jLJlil^  >  ^^®  handled  by  the  servo-multi" 

rx 
plier  units  as  indicated  in  Fig,,  9<. 

A  check  was  made  on  the  analogue  computer  solutions 

obtained  for  the  basic  airplane  (Case  I)  at  sea  level 

using  equations  (111=2^)  of  Ref«  1  for  the  longitudinal 

motion  and  equations  (III-97)  for  the  lateral  motiono 
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